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ABSTRACT 



Context. Intermediate-mass black holes (IMBHs) are of interest in a wide range of astrophysical fields. In particular, the possibil- 
ity of finding them at the centers of globular clusters has recently drawn attention. IMBHs became detectable since the quality of 
observational data sets, particularly those obtained with HST and with high resolution ground based spectrographs, advanced to the 
point where it is possible to measure velocity dispersions at a spatial resolution comparable to the size of the gravitational sphere of 
influence for plausible IMBH masses. 

Aims. We present results from ground based VLT/FLAMES spectroscopy in combination with HST data for the globular cluster NGC 
6388. The aim of this work is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and to compare 
the results with the expected value predicted by the M. - cr scaling relation. 

Methods. The spectroscopic data, containing integral field unit measurements, provide kinematic signatures in the center of the cluster 
while the photometric data give information of the stellar density. Together, these data sets are compared to dynamical models and 
present evidence of an additional compact dark mass at the center: a black hole. 

Results. Using analytical Jeans models in combination with various Monte Carlo simulations to estimate the errors, we derive (with 
68% confidence limits) a best fit black-hole mass of (17 ± 9) x 10 3 Af o and a global mass-to-light ratio of M/L v = (1.6 ± 0.3) M /L g . 

Key words, black hole physics - globular cluster: individual (NGC 6388) - stars: kinematics and dynamics 



1. Introduction 

For a long time, only two mass ranges of black holes were 
known. On the one hand, we have stellar mass black holes, 
which are remnants of massive stars, and can be observed 
in binary systems. On the other hand, there are supermassive 
black holes at the centers of galaxies, some of them accret- 
ing at their Eddington limit and producing the brightest objects 
known (quasars). It has been demonstrated that supermassive 
black holes show a tight correlation between their mass and 
the velocity dispersion of the galaxy in which they reside (e.g. 
Ferr arese & MerritJ l2000t IGebhardt et all l2000t iGultekin et all 
2009). Extrapolating this relation to the lower velocity disper- 
sions of globular clusters, with cr ~ 10-20 kms~' , predicts 
central black holes in these objects with masses of 10 3 - 10 4 M Q . 
Due to the small amount of gas and dust in globular clusters, 
the accretion efficiency of a potential black hole at the center 
is expected to be low. Therefore, the detection of IMBHs at the 
centers of glo bular clusters through X-ray and radio emissions 
is cha llenging (iMiller & Hamiltonll2002r. iMaccarone & Servillaj 
2008). Nevertheless, there is another way to detect IMBHs in 
globular clusters: exploring the kinematics of these systems 
in the central regions. This method, proposed forty years ago 



dBahcall & WolJll976l:IWvUerlll970h , has long been limited by 
the quality of observational datasets, since it requires veloc- 
ity dispersion measurements at a spatial resolution comparable 
to the size of the gravitational sphere of influence for plausi- 
ble IMBH masses (1 - 2" for large Galactic globular clusters). 
However, with existence of the Hubble Space Telescope (HST) 
and with high spatial resolution ground based integral-field spec- 
trographs, the search for IMBHs was revitalized. 

IGebhardt et ail (1 19971 l2000h and iGerssen et al.1 (l2002h 
claimed the detection of a black hole of (3.2 + 2.2) x 1O 3 M 
in the globular cluster M15 from photometric and kinematic 
observations. After more investigations this cluster no longer 



appears as a strong IMBH candidate ( e.g. _ 

Baumgar dTet all 12003. 2005; van den Bos ch et all I2006I) . but 



Dull et al.l fl997l 






* Based on observations collected at the European Organization for 
Astronomical Research in the Southern Hemisphere, Chile (083.D- 
0444). 



new detections of IM BH candidates in other clusters followed. 
Gebh ardt et al.l d2002l 120051) used the velocity dispersion mea- 
sured from integrated light near the center of the M3 1 cluster 
Gl to argue for the presence of a (1.8 + 0.5) x 1O 4 M dark 
mass at the cluster center. The possible presence of an IMBH 
in Gl gained further credence with the dete ction of weak X-ray 
and radio emission from the cluster cen ter dPoolev & Rappapora 
2006: lKongl2007l:IUlvestad et alj|2007t) . Also, the globular clus- 
ter a) Centauri (NGC 5139) has been proposed to host a black 
hole at its center (Novolaet al. 2008, 2010). The authors mea- 
sured the velocity-dispersion profile with an integral field unit 
and used orbit based dynamical models to analyze the data. 
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lAnderson & van der Marell (|2010) studied the same object using 
proper motions from HST images. They found less compelling 
evidence for a central black hole, but more importantly, they 
found a location for the center that differs from previous mea- 
surements. Both Gl and co Centa uri have been suggested to be 
stripp ed nuclei of dwarf galaxies dFreemanll 19931; iMevlan et alj 
2001 ) and therefore may not be the best representatives of glob- 
ular clusters. The key motivation of this work is to probe more 
globular clusters for the presence of IMBHs. 

Further evidence for the existence of IMBHs is the discov- 
ery of ultra luminous X-ray sources at non-nuclear locations 

I II ,J I I 

in sta rburst galaxies (e.g. Fabbiano 1989; Colbert & Mushotzky 
ll999UMatsumoto et al.l l2001; Fabbian oet al.ll200ll) The bright 
est of these compact objects (with L ~ 10 41 ergs -1 ) imply 
masses larger than 1O 3 M assuming accretion at the Eddington 
limit. Several realistic formation scenarios of black holes in 
globular clusters have been developed. The two main forma- 
tion theor ies are: a) IMBHs w ould be Population III stellar 
remnants (Madau & Reea feOOll) . or b) they would form in a 
runaway merging of young stars in sufficiently dense clusters 
(Portegies Zwart et al. 2004; Giirkanetal. 2004; Freitag et al. 
2006). In addition Miller & Hamilton (2002) presented scenar- 
ios for the capture of clusters by their host galaxies and accretion 
in the galactic disk in order to explain the observed bright X-ray 
sources. 

Our goal in this paper is to study the globular cluster NGC 
6388. This cluster is located 1 1.6 kpc away from th e Sun, in 
the outer bulge of our Galaxy. IWhite & Kraf ll {l972) assigned 
it a hi gh metallicity after studying it s color magnitude diagram. 
Later, Illingworth & Freeman (1974) determined the dynamical 
mass of the cluster: with ~ 1.3 x 10 6 M o it belongs to the most 
massive clusters in the Milky Way. In add ition, its high c entral 
velocity dispersion of 18.9 kms 1 dPrvor & Mevlanl fl993), and 
assuming a black-hole mass correlation with velocity dispersion, 
make NGC 6388 a good candidate for detecting an intermediate- 
mass black hole. Besides the kinematic prop erties, the photomet- 
ric cha racteristics are also quite interesting. Noyola & Gebhardt 
(2006, hereafter NG06) found a shallow cusp in the central re- 
gion of the surface brightness profile of NGC 6388. N-body sim- 
ulations showed that this is e xpected for a clu ster hosting an 
intermedi ate-mass black hol e dBaumgardt et al. 2005). For that 
reason, lLanzoni et al.l d2007l) investigated the projected density 
profile and the central surface brightness profile with a combina- 
tion of HST high-resolution and ground-based wide-field obser- 
vations. They found the observed profiles are well reproduced by 
a multimass, isotropic, spherical King model, with an added cen- 
tral black hole with a mass of ~ 5.7 x 1O 3 M . Also, the work of 
Miocchi (2007) suggests the presence of an intermediate-mass 
black hole in NGC 6388 as a possible explanation for the ex- 
tended horizontal branch. 

Another interesting fact about NGC 6388 is t hat it appears to 
conta in multiple stellar populations (e.g. lYoonet al. 2000; Piotto 
2008). This fact makes the scenario of merging clusters plausi- 
ble, leading to potentially complicated dynamics in the center. 
Recently, X-ray observations for NGC 6388 showed no signifi- 
cant signatures for an accreting IMBH dCseh et alj|2010h . But as 
already mentioned, this does not rule out a quiescent black hole. 
In summary, this cluster was chosen as it presents many interest- 
ing features. We measured the kine matics of the central regions, 
allowing us to probe the result of lLanzoni et alj d2007l) with a 
different method, taking the kinematic properties as well as the 
photometric properties into account. 

The basic approach of this work is to first study the light 
distribution of the cluster. Photometric analysis, such as the de- 



Table 1. Properties of the globular clus te r NG C 6388 from 
the re feren ces: NG = Noyola & Gebhardtl d2006l) H= iHarrisI 
( 199 6b. M= iMoretti et alj d2009l) . L= lLanzoni et all d2007l) and 
PM=Prvor & Mevlan ( 1993). 



Parameter 



Value 



Reference 



RA (J200) 


17h 36m 17s 


NG 


DEC (J200) 


-44° 44' 08" 


NG 


Galactic Longitude 1 


345.56 


H 


Galactic Latitude b 


-6.74 


H 


Distance from the Sun ^sun 


11.6 kpc 


M 


Core Radius r c 


7.2" 


L 


Central Concentration c 


1.8 


L 


Heliocentric Radial Velocity V r 


81.2 ± 1.2 km/s 


H 


Central Velocity Dispersion cr 


18.9 km/s 


PM 


Age 


(11.5 ± 1.5) Gyr 


M 


Metallicity [Fe/H] 


-0.6 dex 


H 


Integrated Spectral Type 


G2 


H 


Reddening E(B-V) 


0.38 


M 


Absolute Visual Magnitude M v , 


-9.42 mag 


H 



termination of the cluster center and the measurement of a sur- 
face brightness profile, is described in section [2] De-projecting 
this profile gives an estimation of the gravitational potential pro- 
duced by the visible mass. The next step is to study the dynamics 
of the cluster. Section|3]gives an overlook of the FLAMES obser- 
vations and data reduction and section [4] describes the analysis 
of the spectroscopic data. With the resulting velocity-dispersion 
profile, it is possible to estimate the actual dynamical mass. The 
next step is to compare the data to dynamical Jeans models (sec- 
tion[5]l. These models take the light profile and predict a velocity- 
dispersion profile, which is scaled to the data in order to obtain 
the mass-to-light ratio. This is done for models containing dif- 
ferent black-hole masses until the best fit to the observed profile 
is found. In section [6] we summarize our results, list our conclu- 
sions and give an outlook for further studies. 



2. Photometry 

The photometric data, retrieved from the archive, were ob- 
tained with the Advanced Camera of Surveys (ACS) of the 
Hubble Space Telescope (HST) in the Wide-Field Channel (un- 
der the HST program SNAP-9821, PI: B.J. Pritzl) between 
October 2003 and June 2004. This data set is composed of six 
B (F435W), V (F555W) and I (F814W) images with exposure 
times of 11, 7 and 3s, respectively. It gives a complete cover- 
age of the central region of the cluster out to a radius of 110". 
The data were calibrated, geometrically corrected and dither- 
combined as retrieved from the European HST-Archive (ST- 
ECF, Space Telescope European Coordinating Facilityfl). 

2.1. Color magnitude diagram (CMD) of NGC 6388 

The CMD was obtained using the programs daophot II, allstars 
and airframes by P. Stetson, applied to our HST image. F or a de- 
tailed documentation of these routines see lStetsonldl987l) . These 



1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science 
Institute (STScI/NASA), the Space Telescope European Coordinating 
Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 
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programs were especially developed for photometry in crowded 
fields and therefore ideally suited for the analysis of globular 
clusters. The routines find, phot and psf identify the stars, per- 
form aperture photometry and compute the average point spread 
function (PSF) over the image, respectively. 

Once the PSF has been defined, the next step is to group the 
neighboring stars to apply the multiple-profile-fitting routine si- 
multaneously by the task allstars. Afterwards, the find task is 
applied again to find, in the star-subtracted image, stars which 
were not found in the first run. The entire procedure was per- 
formed on the V- and I-band images independently. As a next 
step the programs daomaster and allfram.es were used to com- 
bine both images and to create the final catalog containing all 
stars, their positions and magnitudes in the two bands. At the 
end, we calibrated the final catalog to the Johnson magnitude 
system by following the steps described in Sirianni et al. (2005). 

In order to get better quality at the faint end of the CMD, 
one final step was applied to the catalog. The program separa- 
tion (iBruntt et al.ll20 03) computes a separation index for every 
star in a catalog. This index is calculated by, first, evaluating a 
local surface brightness at the position of a given star. Second, 
the local surface brightness is compared to the sum of the surface 
brightnesses produced by the PSF of all the other stars in the field 
at the position of the centroid of that star. The ratio of these two 
surface brightness values expressed in magnitudes determines 
the separation index. Thus, the stars could be selected consider- 
ing background-light contamination and not only by magnitude. 
Figure Q] shows the final CMD of stars with a separation index 
> 5 overplotted with the positions of the brightest stars in the 
ARGUS pointing and the spectroscopic template star (see sec- 
tion|3). 

2.2. Center of the cluster 

Using the star catalog generated with daophot, the center of 
the cluster can be determined. Precise knowledge of the clus- 
ter center is important since the shape of the surface bright- 
ness and the angular averaged line-of-sight velocity distribu- 
tion (LOSVD) profiles depend on the position of that cen- 
ter. Using the wrong center typically produces a shallower in- 
ner profile. For example determining the c enter for a> Centauri 
is not trivial (lNovola& Gebhardt 200^ iNovola et all l2010t 
lAnderson &va n der Marel 2010). Fortunately, NGC 6388 is not 
as extended and has a steeper light profile than a> Centauri so 
that the center is easier to determine. NG06 determined the cen- 
ter of this cluster to be at a - 17 : 36 : 17.18, 5 - -44 : 
44 : 07.83 (72000), with an uncertainty of 0.5", by minimizing 
the standard deviation of star counts in eight segments of a cir- 
cle. In view of the discrepancies about the center location for u 
Centauri, we decided to recompute the center of NGC 6388 and 
to evaluate how precisely the center can be determined. NG06's 
center on our I-band image was used as a first guess for the fol- 
lowing methods to determine the center in the ACS images, i.e. 
our reference frame for the spectroscopy. 

The first method is a simple star count as described in 
McLaughlin et al. (2006). The catalog generated with daophot 
(see section l2~TT i contains 88,406 stars. In a field of 4" x4" a grid 
of trial centers was created, using a grid spacing of 2 ACS pixels 
(0.1"). Around each trial center a circle of 300 pixels (15") was 
considered and divided into 16 wedges as shown in Figure|2] The 
stars in each wedge were counted and compared to the opposite 
wedge. The differences in the total number of stars between two 
opposing wedges were summed for all 8 wedge pairs. The coor- 
dinates that minimized the difference defined a first guess of the 
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Fig. 1. The color-magnitude diagram of NGC 6388. Overplotted 
are the brightest stars identified in the ARGUS field of view (red 
circles), and the template star used (star symbol). 



center of the cluster. This center was refined as described in the 
next sections. We compared our result to the center obtained by 
NG06. The two centers are only 0.32" apart and thus coincide 
within the error bars of 0.5" (as determined by NG06 performing 
artificial image tests). 

The second method that we used is also described in 
iMcLau^hlin et al.l d2006l) . Instead of comparing the total num- 
ber in pairs of opposite wedges, a cumulative radial distribution 
for each wedge in 4 bins was generated. This time 8 wedges were 
used instead of 16 to avoid too large stochastic errors due to in- 
sufficient numbers of stars in each bin. The bins were placed 
at equidistant radii. Again the absolute value of the integrated 
difference between the radial distributions in any two opposing 
wedges was calculated and the minimum used to determine the 
new center of the cluster. The so derived center lies within 0.1" 
of the one derived with the previous method and within 0.3" of 
NG06's location. 

We present a last method in which the light of the stars in- 
stead of their number, is considered. Similar to the first method, 
16 wedges without any radial bins were generated, but this time 
not the stars were counted but the luminosities of the stars were 
summed up in each wedge. In order to avoid a bias by the con- 
tamination of a few bright stars, only stars from the lower gi- 
ant branch and the horizontal branch (between my = 15 and 
my - 19) were used. This approach reproduces NG06's center 
within the error bars (0.28") as well. The method is illustrated in 
Figure [2] Shown are the grid of trial centers, the 16 wedges we 
applied to each of them, and the final contour plot. 
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Fig. 2. The method to determine the globular cluster center. From left to right: a zoom into our HST/ACS image of NGC 6388 
(used to derive the CMD and luminosity profile for the cluster) with the over-layed grid of trial centers. Around each trial center, 
differential stellar counts and stellar luminosity were computed in two opposite wedges extending 15" in radius [middle panel]. A 
minimization of the residuals, as shown in the contour plot [right panel] determines the center to within 0.2". The distances between 
the contours vary between 40 - 70 stars for the star count methods. 



For the subsequent analysis, we used the center derived from 
our ACS catalog. To estimate the error of our center determina- 
tion we took the scatter from the three different methods as well 
as an additional test where only eight wedges where used. For 
these eight wedges we repeated the routine by only using the 
cardinal wedges and the semi-cardinal wedges separately. For 
this dataset (HST J8ON08OYQ, also used as position reference 
coordinates) we derived a final position of the center of: 



(x c ,y c ) = (2075.5,2548.5) + (3. 1,2.1) pixel 



a = 17 : 36 : 17.441, Aa = 0.2" (J2000) 
S = -44:43 : 57.33, A6 = 0.1" 



(1) 

(2) 
(3) 



All the derived centers lie within a few tenths of an arcsecond 
(~ 10~ 2 pc) radius. However, it must be considered that all the 
methods are using the same catalog derived by DAOPHOT This 
catalog most likely suffers from incompleteness since the bright 
stars are covering the fainter stars underneath. This could bias 
the center towards the bright stars even if they are excluded in 
the counts. Our center agrees within the error bars with the one 
determined in NG06 on WFPC2 images. It also coincides with 
the center found by lLanzoni et al.l (l2007t> which, according to 
them, lies ~ 0.5" northwest of the NG06 center and therefor 
closer to ours. 



2.3. Surface brightness profile 

The last step of the photometric analysis was to obtain the sur- 
face brightness profile. This is required as an input for the Jeans 
models described in the following section. To obtain the profile, 
a simple method of star counts in combination with an integrated 
light measurement from the ACS image was applied. The fluxes 
of all stars brighter than my = 18 were summed in radial bins of 
equal width (50 pix) around the center and divided by the area 
of each bin. Since there are no stars in the gap between the two 
ACS CCD chips this area was subtracted from each affected bin. 



In addition, the integrated light for stars fainter than my = 18 
was measured directly from the HST image. Using the same ra- 
dial bins as in the star count method, we measured the statistical 
distribution of counts per pixel excluding regions with stars with 
my < 18. After trying different methods to derive the "average" 
of the distribution of counts, we applied a simple average that 
takes into account the faintest pixels i.e. stars. At the end, the 
flux per pixel was transformed back into magnitude per square 
arcseconds and added to the star counts profile. 



L, 



We compared 



anzoni 



compared 

etal. (2007) 



our profile with Trag er et alj d 1995b . 

and NG06. The latter was derived by mea- 
suring the integrate d light using a bi-weight estimator while 
lLanzoni et al.l d2007l) derived their points by taking the average 
of the counts per pixel in each bin. We were able to reproduce the 
shape in the outer regions, but due to the method and data that we 
used our profile in the innermost region has a high uncertainty. 
The errors on our profile were obtained by Poisson statistics of 
the number of stars in each bin. With a linear fit inside the core 
radius (~ 7") we derived a logarithmic slope of /3 — 0.7 + 0.2 
(where /3 corresponds to [iy ~ log r 13 with the surface brightness 
jUy), which results in a slope of the surface luminosity density 
I(r) a r a of a = -0.28 ± 0.08. This value is steeper (but con- 
sistent within the errors) than the slope of a — -0.13 + 0.07 
derived by NG06 and also consistent with the slope derived by 
lLanzoni et all d2007l) of a = -0.23 ± 0.02. 

To derive the SB profile, we used a calibrated catalog but we 
also combined this with measurements that we obtained from 
counts per pixel. It is possible that the transformations from 
counts into magnitudes do not match the calibrations of the cat- 
alog and therefore would cause an offset in our profile. The fact 
that we used all the light in the image, including the background 
light could also shift the profile to higher magnitudes than the 
intrinsic brightness. For that reason, we scaled o ur profile to 
the outer parts of the profile bv lTragere t al. ( 1995). This profile 
was obtained with photometrically calibrated data and provides 
a good reference for the absolute values of the profile. The final 
result of the surface brightness profile is shown in Figure [3] 
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Table 2. The derived surface brightness profile in the V-band. 
AV/, and AV/ are the high and low values of the errors, respec- 
tively. 



logr 
[arcsec] 



V 



AV, 



AV,, 



[mag/arcsec 2 ] [mag/arcsec 2 ] [mag/arcsec 2 ] 



0.40 


14.00 


0.52 


0.35 


0.12 


14.53 


0.52 


0.35 


0.18 


14.31 


0.21 


0.18 


0.40 


14.74 


0.16 


0.14 


0.60 


14.78 


0.10 


0.09 


0.81 


15.18 


0.07 


0.07 


1.06 


15.62 


0.05 


0.05 


1.22 


16.21 


0.05 


0.05 


1.33 


16.68 


0.06 


0.06 


1.42 


17.18 


0.07 


0.07 


1.50 


17.50 


0.07 


0.07 


1.56 


17.91 


0.08 


0.07 


1.62 


18.20 


0.08 


0.08 


1.67 


18.54 


0.09 


0.09 


1.71 


18.74 


0.10 


0.09 


1.75 


18.80 


0.09 


0.09 


1.79 


19.17 


0.11 


0.10 


1.82 


19.33 


0.11 


0.10 


1.85 


19.61 


0.13 


0.11 


1.88 
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0.10 
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Fig. 4. The positions of the three ARGUS pointings (A, B and 
C) reconstructed on the HST/ACS image. 



three different positions (exposure times of the pointings: A: 
3 x 480s + 3 x 1500s, B: 3 x 1500s, C: 1 x 1500s) to cover 
the sphere of influence of the potential black hole (see Figure|4]i. 
The IFU fibres were placed around the core radius of the cluster 
to obtain velocity dispersions in the outer regions. 

The kinematics were obtained from the analysis of the 
Calcium Triplet (~ 850 nm) which is a strong feature in the 
spectra. The expected velocity dispersions lie in the range 5- 
20 kms -1 and need to be measured with an accuracy of 1-2 
kms 1 . This implied using a spectral resolution around 10000, 
available in the low spectral resolution mode set-up LR8 (820 - 
940 nm, R= 10400). 




0.5 1.0 

log r (arcsec) 

Fig. 3. The surface brightness profile of NGC 6388. The profile 
shows a clear cusp for the inner 10 arcseconds. Also shown is the 
MGE fit (solid line) which was used to parametrize our profile 
(see section[8]) 



3. Spectroscopy 

3.1. Observations 

The spectroscopic data were observed with the GIRAFFE spec- 
trograph of the FLAMES (Fiber Large Array Multi Element 
Spectrograph) instrument at the Very Large Telescope (VLT) 
in ARGUS (Large Integral Field Unit) and IFU (Integral Field 
Unit) mode. The set contains spectra from the center (ARGUS) 
and the outer regions (IFU) for the globular cluster NGC 6388. 
The observations were performed during two nights (2009-06- 
14/15). The ARGUS unit was set to the 1 : 1 magnification 
scale (pixel size: 0.3", 14 x 22 pixel array) and pointed to 



3.2. Data reduction 

We reduced the spectroscopic data with the GIRAFFE pipeline 
programmed by the European Southern Observatory (ESO). 
This pipeline consists of five recipes which are briefly described 
below. 

First, a master bias frame was created by the recipe gimas- 
terbias from a set of raw bias frames. Next, a master dark was 
produced by the recipe gimasterdark which corrected each in- 
put dark frame for the bias and scaled it to an exposure time 
of 1 second. The recipe gimasterflat was responsible for the de- 
tection of the spectra on a fiber flat-field image for a given fiber 
setup. Gimasterflat located the fibers, determined the parameters 
of the fiber profile by fitting an analytical model of this profile to 
the flat-field data and created an extracted flat-field image. This 
image was later used to apply corrections for pixel-to-pixel vari- 
ations and fiber-to-fiber transmission. 

The pipeline recipe giwavecalibration computed a disper- 
sion solution and a slit geometry table for the fiber setup in use. 
This was done by extracting the spectra from the bias corrected 
ThAr lamp frame, using the fiber localization (obtained through 
the flat-field), selecting the calibration lines from the line cata- 
log, and predicting the positions of the ThAr lines on the detec- 
tor using an existing dispersion solution as an initial guess. The 
quality of the calibration was checked in two ways. The first one 
was by measuring the centroid of one dominating sky emission 
line in all 14 sky spectra and the second was a cross-correlation 
in Fourier-space of all sky spectra. 

Both methods did not show any systematic shift of the spec- 
tra and resulted in an RMS of 0.03 A, which corresponds to a 
velocity of 1 km s -1 . The last recipe, giscience combined all cal- 
ibrations and extracted the final science spectra. A simple sum 
along the slit was applied as extraction method. The input obser- 
vations were averaged and created a reduced science frame, the 
extracted and rebinned spectra frame. At the end, the recipe also 
produced a reconstructed image of the respective field of view 
of the IFU and the ARGUS observations. 
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4. Kinematics 

This section describes how we created a velocity map in order 
to check for rotation or other peculiar kinematic signals and how 
we derived a velocity-dispersion profile. This profile is then used 
to fit analytic models, described in the next section. 



4.1. Velocity map 

The reconstructed images of the three ARGUS pointings were 
matched to the HST image. With this information, the pointings 
were stitched together and each spectrum correlated to one posi- 
tion in the field of view. The resulting catalog of spectra and their 
coordinates allowed us to combine spectra in different bins. The 
combined pointing contains 24 x 29 spaxels. 

For each spaxel, the penalized pixel-fitting (pPXF) program 
developed by Cappellari & Emsellem (2004) was used to derive 
the kinematics in that region. Figure [6] shows a) the field of view 
of the three pointings on the HST image, b) the field convolved 
with a Gaussian and resampled to the resolution of the ARGUS 
array and c) the reconstructed and combined ARGUS image. 
The corresponding velocity map is shown in Figure [7] For ev- 
ery ARGUS pointing, a velocity map was derived before com- 
bining. The C pointing is a bit more noisy due to the fact that 
it only had one exposure. To check for systematic wavelength 
offsets we compared the derived velocities of the different point- 
ings and exposures at overlapping spaxels. We found offsets of 
up to 1 km s. We checked for systematic wavelength offsets 
between the pointings by cross-correlating their averaged sky 
spectra. The rms of the shifts measures 0.006 pixel (0.0018 A) 
which corresponds to a velocity shift of ~ 0.1 kms.~' We con- 
clude that the velocity offsets are not due to uncertainties in the 
wavelength calibration and do not cause problems for our further 
analysis. Instead, the offsets are explained as follows: 

i) Different pointings were dithered about half an ARGUS 
spaxel. This means that in every pointing a slightly different 
combination of stars contributes to the different spaxel, which 
could cause a different velocity measurement, ii) Errors in the 
reconstruction of the pointing positions could also cause veloc- 
ity offsets between two positions. If the spaxels do not point at 
the exact same position we measure, similar to point i), differ- 
ent velocities, iii) Not all pointings had the same exposure time. 
This causes different signal-to-noise ratios for different point- 
ings. This could also explain some variations in the velocity 
measurements with the pPXF routine. Note that for the veloc- 
ity dispersion measurements, we grouped the individual spectra 
in large annuli. For that reason, the position offsets do not cause 
a problem for the further analysis as the "error" in position is 
not propagated into the velocity dispersion measurement as a 
systematic shift. 

We used a star from the central pointing (pointing A) as a 
velocity template. This has the advantage that it went through 
the same instrumental set-up as well as the same reduction steps 
as all other spectra and therefore the same sampling and wave- 
length calibration. The green circle in the HST image in Figure 
[6] (left panel) and the blue star in the color magnitude diagram 
(Figure[TJ marks the star which was used. We also identified the 
brightest stars form the pointing in the CMD to make sure that 
the template and other dominating stars are not foreground stars 
(see Figure [TJ. In order to derive an absolute velocity scale, the 
line shifts of the templates were measured by fitting a Gaussian 
to each line and deriving the centroid. This was compared with 
the values of the Calcium Triplet in a rest frame and the average 



Fig. 5. Combined spectra of the first, third and sixth bin over- 
plotted by their best fit (red line). Due to the binning, the spectra 
show a high signal-to-noise ratio. 



The most important parts of our reduction are the sky sub- 
traction and an accurate wavelength calibration to avoid arti- 
ficial line broadenings due to incorrect line subtractions. The 
program we used was developed by Mike Irwin and described 
in Battagl ia et alj Q2008). To subtract the sky, the program first 
combined all (14) sky fibers using a 3-sigma clipping algorithm 
and computed an average sky spectrum. Using a combination 
of median and boxcar, it then split the continuum and the sky- 
line components in order to create a sky-line template mask. For 
the object spectra, the same method of splitting the continuum 
from the lines was applied. In the process, the sky-lines were 
masked out to get a more accurate definition of the continuum. 
Afterwards, the sky-line mask and the line-only object spectra 
were compared to find the optimum scale factor for the sky spec- 
trum in order to subtract the sky-lines from the object spectra. 
The continuum was added back to the object spectra and the sky 
continuum subtracted by the same scaling factor as obtained for 
the lines (assuming that lines and continuum have the same scal- 
ing factor). 

After applying the sky subtraction program, the last step in 
the reduction of the data was to remove the cosmic rays from the 
sp ectra. This wa s done using the program LA-Cosmic developed 
by van Dokkum ( 2001 ) and based on a Laplacian edge detection. 
In order to avoid bright stars dominating the averaged spectra 
when they were combined, we applied a simple normalization 
by fitting a spline to the continuum and dividing by it. 

For the small IFUs, we applied the same reduction steps 
as for the large integral-field unit, except for the cosmic ray 
removal, since LA-Cosmic did not give a satisfying result. 
However, we were able to average all exposures applying a 
sigma clipping method in order to remove the cosmic rays, due 
to the fact that these pointings were not dithered. 
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Fig. 6. The ARGUS field of view: the HST image is shown on the left. The red circle marks the center, the green one the template 
star which was used. The same image but convolved with a Gaussian and resampled to the 0.3 pixel scale of the ARGUS array 
[center]. Compared to the actual reconstructed ARGUS pointing [right], it is clear that they are pointing to the same region. 
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Fig. 7. The velocity map of NGC 6388. As shown on the velocity scale, the blue spaxels indicate approaching velocities and the red 
ones receding. The white stars mark the spaxels which we excluded from the velocity dispersion measurement as they might suffer 
from shot noise. Also shown are the first three velocity dispersion bins to help visualize the binning method. 



shift was calculated. We then corrected the radial velocity for the 
heliocentric reference frame. 



4.2. Shot noise corrections 

An obvious feature of the velocity map in Figure Q is the bright 
blue spot in the middle left of the pointing which seems to dom- 
inate the blue (approaching) part of the map. We investigated 
whether we see a real rotation around the center or just one or a 
few bright stars with a peculiar velocity dominating their envi- 



ronment and with their light contaminating neighboring spaxels. 
Therefore, we tested the influence of each star on the adjacent 
spaxels. This allowed us to test if the derived velocity disper- 
sion is representative of the entire population at a given radius or 
whether it is biased by a low number of stars, i.e. dominated by 
shot noise. 

To perform this test, we considered our photometric cata- 
log (described in section 13 for the field of view covered by 
the ARGUS pointings. At every position of a star in the cata- 
log, a two dimensional Gaussian was modeled with a standard 
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Table 3. The four moments of the velocity distribution of NGC 
6388. 



logr 


V 


Vrms 


Hi 


M 


[arcsec] 


[km/s] 


[km/s] 






-0.05 


-1.6 ±2.9 


23.3 ±3.1 


-0.05 ±0.16 


-0.16 ±0.13 


0.26 


-2.2 + 0.6 


25.7 ±1.8 


-0.07 ± 0.01 


-0.10 ±0.01 


0.43 


-0.1 ±0.6 


24.0 ±1.3 


-0.08 ± 0.02 


-0.07 ± 0.02 


0.56 


2.9 ±0.6 


22.2 ±1.1 


-0.09 ± 0.02 


-0.09 ± 0.02 


0.65 


6.4 ±0.6 


20.5 ± 0.9 


-0.08 ± 0.02 


-0.07 ± 0.02 


0.78 


4.6 ±0.8 


18.6 ±0.8 


-0.03 ± 0.04 


-0.01 ± 0.03 


1.04 


3.5 ±0.8 


16.8 ±2.9 


-0.12 ±0.03 


-0.14 ±0.02 


1.26 


6.7 ±1.0 


13.6 ±2.9 


-0.07 ± 0.04 


-0.03 ± 0.04 



deviation set to the seeing of the ground based observations 
(FWHM = 0.9") and scaled to the total flux of the star. The 
next step was to measure the absolute amount and fraction of 
light that each star contributes to the surrounding spaxels. After 
computing these values for every star in the pointing, we had 
the following information for each spaxel: a) how many stars 
contribute to the light of that spaxel, and b) which fraction of the 
total light is contributed by each star, i.e. we determined whether 
the spectrum in a given spaxel was dominated by one or a few 
stars. The test showed that most of the spaxels contained mean- 
ingful contributions by more than 10 stars. Some spaxels, how- 
ever, were dominated by a single star contributing more than 50 
% to the spaxel's light. For this reason, the contribution in per- 
cent of the brightest star was also derived by the program. We 
found out that the blue area in the left side of the velocity map of 
Figure|7]is not due to a single star. In fact this area in the velocity 
map corresponds to a group of at least 10 stars moving with 10 - 
40 kms -1 with respect to the cluster systemic velocity. 

4.3. Velocity-dispersion profile 

To derive a radial velocity-dispersion profile for the stellar pop- 
ulation of NGC 6388 it is necessary to bin the spectra accord- 
ingly. We divided the pointing into six independent angular bins, 
each of them with the same width of three ARGUS spaxels 
(0.9", 0.04 pc). To check the effect of the distribution of the bins 
on the final result, we tried different combinations of bins and bin 
distances as well as overlapping bins. We found no change in the 
global shape of the profile when using different bins. Therefore 
the choice of the bins was not critical, but in order to make an 
accurate error estimation, independent bins are more useful. 

In each bin, all spectra of all exposures where combined with 
a sigma clipping algorithm to remove any remaining cosmic 
rays. Velocity and velocity-dispersion profiles were computed 
using the pPXF method applied to the binned spectra. The ve- 
locity map in Figure [7] shows the dynamics in the cluster center. 
A rotation-like or a shearing signature is visible, but it seems 
to be a very local phenomenon (within 3", 0.15 pc). In earlier 
attempts, we split the pointing in two halves (tilted line in fig- 
ure 13 in order check for consistency and symmetry on both 
sides. The velocity dispersion was then derived separately for 
the binned spectra on each side. Both sides show a rise in the 
velocity-dispersion profile but differ in their shape and absolute 
values from each other. However, treating both sides separately 
would not properly take into account a possible rotation at the 
center and therefore not exactly represent the observed data. We 
decided to use the total radial profile over the 360 degree bins 

and measure the second moment Vrms = -v/V^, + cr 2 , with the 



rotation velocity V rot and the velocity dispersion cr. The reason 
why we chose to analyze Vrms instead of cr is twofold. First, the 
Jeans models require the Vrms rather than the pure velocity dis- 
persion as an input. Second, the broadening of the line we mea- 
sured represents Vrms and not cr. The velocity dispersion can be 
obtained by subtracting the rotation velocity from this quantity. 
Determining the rotation is difficult due to the large shot noise 
introduced by the small number of spatial elements in the cen- 
tral region. The second moment, however, is robustly measured 
since we average over all angles. For simplicity we refer to the 
Vrms profile as the velocity-dispersion profile in our study. 

In addition to the central pointings, we derived kinematics 
for regions further out using the small IFU measurements, which 
were scattered at larger radii around the cluster. Unfortunately, 
the surface brightness of the cluster drops quickly with radius 
resulting in a low signal-to-noise ratio for the IFUs most distant 
from the center. These could therefore not be used for further 
analysis. Only the two innermost positions showed reasonable 
signal, so that these two pointings could be used as two separated 
data points at 1 1" and 18" radius, respectively. The disadvantage 
of the small IFU fields (20 spaxels for a total field-of-view of 
3" x 2") is that these values are very affected by shot noise since 
only a few stars fall into the small field-of-view. Consequently, 
these points show larger errors than the rest of the profile. 

Further, we estimated the radial velocity of the cluster in 
the heliocentric reference frame. We combined all spectra in the 
pointing and measured the velocity relative to the velocity of the 
template. This value was then corrected for the motion of the 
template and the heliocentric velocity and results in a value of 
V r = (80. 6 + 0.5) km s' 1 which is in good agreement with the 
value from lHarrisldl996l) V r = (81.2 ± 1 .2) km s _I . 

We ran Monte Carlo simulations to estimate the error on 
Vrms due to shot noise. From the routine described in section 
14.21 we knew exactly how many stars are contributing (nStars) 
and how many spaxels are summed (nSpax) in each bin. We took 
all stars detected in each bin and their corresponding magnitude. 
Each of the stars was then assigned a velocity chosen from a 
Gaussian velocity distribution with a fixed dispersion. We used 
our template spectrum and created nSpax spaxel by averaging 
nStars flux weighted and velocity shifted spectra. The result- 
ing spaxels were then normalized, combined and the kinematics 
measured with pPXF (as for the original data). After 1000 real- 
izations for each bin, we obtained the shot noise errors from the 
spread of the measured velocity dispersions. For the outer IFU 
pointings, we extrapolated the surface density to larger radii and 
performed the same Monte Carlo simulations as for the inner 
pointing, with random magnitudes drawn from the magnitude 
distribution in this region. The errors for the velocity and the 
higher moments were derived by applying Monte Carlo simula- 
tions to the spectrum itself. This was done by repeating the mea- 
surement for 100 differen t realizations of data by adding noise to 
the original spectra, (see Cappellari & Emsellem 2004, section 
3.4) 

The resulting profile is displayed in Figure[8] The plot shows 
the second moment of the velocity distribution Vrms- Except for 
the innermost point, a clear rise of the profile towards the center 
is visible. The highest point reaches more than 25 km s before 
it drops quickly below 20 km s at larger radii. In table[3]we list 
the results of the kinematic measurements. The first column lists 
the radii of the bins. The following columns show the central 
velocities of each bin in the reference frame of the cluster, the 
second moment Vrms as well as higher velocity moments h3 and 
h4. It is conspicuous that all h4 moments tend to have negative 
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values. This hints at a lack of radial orbits in the central region 
of the globular cluster. 

For general purposes we also determined the effectiv e veloc- 
ity dispersion cr e as described in iGtiltekin et al.l d2009t) . This is 
defined by: 



<r? 






VUKr)dr 



f Q R 'I(r)dr 



(4) 



Where R e stands for the half light radius (~ 40") and I(r) the 
surface brightness profile. We extrapolated our kinematic data to 
the half light radius since our furthest out data point only reaches 
a radius of 18". This results in an effective velocity dispersion of 
cr e - (18.9 + 0.3) km s" 1 which is in per fect agreement with the 
value derived bv lPrvor & Mevlanl ([993) listed in tableQ] 



the contour plot ofy 2 as a function of black-hole mass and mass- 
to-light ratio over a grid of isotropic models (black points). The 
contours represent Ay 2 = 1 .0, 2.7, 4.0 and 6.6 which correspond 
to a confidence of 68, 90, 95 and 99 percent for 1 degree of free- 
dom (marginalized). This implies that for an isotropic model and 
this specific surface brightness profile the models predict a black 
hole of at least M, — 5 x 1O 3 M with a confidence of 99 percent. 
We also tested whether anisotropic Jeans models would re- 
sult in a significantly better fit. To do this, we repeated the model 
fitting with j3 + values. We found a decrease of x 2 for a ris- 
ing /? down to^- 2 = 4.21 for/? = 0.5 with a lower black-hole 
mass of M, - 5.3 x 1O 3 M . However, this requires a very high 

anisotropy of /3 — 1 - Vg/V 2 = 0.5 which is expected to be un- 
stable in the center of a relaxed globular cluster such as NGC 
6388. The anisotropy will be discussed in more detail in section 



5. Dynamical models 

The main goal of this work is to compare the derived kinematics 
and light profile with a set of simple analytical models in order 
to test whether NGC 6388 is likely to host an intermediate-mass 
black hole in its center. We used J eans models as implemented 
and described in Cappellaril (120081) . 



5.1. Isotropic and anisotropic models 

The first input for the Jeans models is the surface brightness pro- 
file in order to estimate the 3-D density distribution in the clus- 
ter. Given the fact that the density can only be observed in a pro- 
jected way, the profiles have to be deprojected. One way of doing 
this is by applying the multi-Gaussian expansion (MGE) method 
developed by Emsellem et al. ( 1994). The basic approach of this 
method is to parametrize the projected surface brightness with a 
sum of Gaussians since the deprojection of a Gaussian function 
results again in a Gaussian. 

To apply this parametrization and to compare results of the 
Jeans equation to our data, we used the Jeans Anisotropic MGE 
(JAM) dynamical models for stellar kinematics of spherical and 
axisymmetric galaxies, as well as the multi-Gaussian expansion 
implementation developed by Cappellari (2002, 2008). The IDL 
routines provided by M. Cappellaro enable the modeling of the 
surface brightness profile and fitting the observed velocity data 
and the mass-to-light ratio at the same time. We used our surface 
brightness profile in combination with a spherical model with 
different values of anisotropy and constant M/Ly setups along 
the radius of the cluster. 

In Figure[8] we plot the isotropic model data comparison for 
our velocity-dispersion profile. The JAM program does not ac- 
tually fit the model to the data but rather calculates the shape of 
the second moment curve from the surface brightness profile, de- 
convolves the profile with the PSF of the IFU observations, and 
then scales it to an average value of the kinematic data in order 
to derive the MjLy. This explains why all trial models meet in 
one point. The thick black line marks the model with the lowest 
X 2 value and therefore the best fit. 

With the given surface brightness profile and without a cen- 
tral black hole, the models predict a drop in velocity disper- 
sion towards the central region. As a final result, we used the 
X 2 statistics of the fit to estimate an error. This results in: M. - 
(18 ±6) xl0 3 M o and M/L v = (1.7 + 0.2) M Q /L G . Figure|9]shows 



5.2. Error estimation 



In section l231 we introdu ced the different sur face brightness pro- 
files for NGC 6388. The iTrager etail d 1 995l) profile contains the 
most data points and extends to large radii (~ 4')- However, the 
inner regions (important for our dynamical purposes) are not as 
well sampled by the ground based observations used by Trager 
et al. The profile by NG06 was derived by measuring integrated 
light on a WFPC2 image using a bi-weight estimator and covers 
the inner regions of the cluster very well. To calibrate the profiles 
to the correct absolute values NG06 also adjusted it to the Trager 
et al. profile. The third profile was obtained by lLanzoni et al.l 
(2007) by computing the average of the counts per pixel in each 
bin of their ACS/HRC images. Since it was not clear in the pa- 
per how they calibrated their photometry, this profile was also 
shifted (by a small amount of 0.09 mag) to overlap with the 
Trager et al. profile. As a test, we used each of these profiles 
as an input for the Jeans models and found similar black-hole 
masse s varying from M, - (25 + 6) x 1O 3 M for the Trager et al 
(1 19951) profile to M. - (11 ± 5) X 1O 3 M for the Lanzo ni et al 



Available at http://purl.org/cappellari/idl 



(2007) data. The result from the profile provided by NG06 is 
identical to the fit of our profile in black-hole mass and slightly 
lower in the mass-to-light ratio (M/L v - (1.5 ± 0.2) M Q /L Q ). 
Again we find the lowest x 1 value (and therefore the best fit) 
for an anisotropic model of/? = 0.5. In this case the observations 
can be explained without any black hole for all profiles. As men- 
tioned in the previous section this rather unlikely configuration 
will be discussed in more detail in the next subsection. 

From our studies of the surface brightness profiles, we see 
that they yield similar, but not identical results. The shape of 
the SB profile predicts the shape of the velocity-dispersion pro- 
file. Therefore, it is crucial to test the effect of variation of the 
surface brightness profile in the inner regions. To perform this, 
we run Monte Carlo simulations on the six innermost points 
of the profile. We used 1000 runs and a range of/? between 
-0.5 and 0.5. The results are displayed in figure [TOl The black 
points mark the mean value of the derived black-hole masses. 
The shaded contours represent the 68,90 and 95 percent confi- 
dence limits. The mean black-hole mass decreases at higher ft 
values since a radial anisotropy can mimic the effect of a cen- 
tral black hole. Nevertheless, for an isotropic case, we have a 
black hole detection within a mass range of 17^ x 1O 3 M and 
a global mass-to-light ratio of MjLy - 1.6*9? M e /L Q (errors 
are the 68 percent confidence limits). We derived the total error 
(resulting from kinematic and photometric data) by again, run- 
ning Monte Carlo simulations and varying both, the velocity- 
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Fig. 8. Different isotropic (J3 - fixed) Jeans models applied to the kinematic data from our pointing. The M/L is fitted to the data 
by scaling the profile to the data points. The final values are obtained by finding the minimum of x 2 while varying the black-hole 
mass: j3 - 0.0, M BH = 18.24 x 1O 3 M , M/L v - 1.65, xl lin = 8.89. The plot on the left panel shows the models together with the 
data and the best fit. Labeled are also the black holes masses of the two enclosing models. The solid black line marks the best fit. 
The right panel shows the^ 2 values of every model as well as the Ay 2 = 1 line (dashed line). 
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Fig. 9. The contours of the x 2 as a function of black-hole 
mass and mass-to-light ratio. Each point represents a particu- 
lar isotropic model. The plotted contours are Ax 2 - 1.0,2.7,4.0 
and 6.6 implying a confidence of 68, 90, 95 and 99 percent. 



dispersion profile and the surface brightness profile at the same 
time. As a final result for the uncertainties of our method we de- 
rived SM. = 9x 1O 3 M , 6M/L V = 0.3 M o /L . This shows, 
that a significant percentage of the error (~ 75%) in the context 
of isotropic Jeans models results from the uncertainties of the 
surface brightness profile. 

5.3. Anisotropy in relaxed globular clusters 

In the previous sections we have shown, that the strong signif- 
icance for a black hole with isotropic models vanishes when 
using anisotropic models. The anisotropic Jeans models show 
a better fit for a model with much lower black-hole masses or 
without a black hole requiring an anisotropy of j3 - 0.5 or higher. 
Thus, it is necessary to investigate the probability of such a kine- 



matic configuration in our cluster. We run N-body simulations 
for 30000 equal-mass stars, distributed according to a Plummer 
model. The initial anisotropy was created by a cold collapse, 
which resulted in a cluster that was slightly anisotropic in its core 
and had increasing anisotropy further out rising to /3 - 0.9 out- 
side the half-mass radius. Furthermore, we took the tidal field of 
the galaxy into account by letting the cluster move on a circular 
orbit around a point-mass galaxy. The ratio between tidal radius 
and core radius was set to 1:10 as it is the case for most globular 
clusters as well as for NGC 6388. We then let the cluster evolve 
over a time scale of several relaxations times and calculated the 
anisotropy for snapshots of the model. 

Figure [TT] shows the result of these tests. The different col- 
ors represent different areas in the cluster. Important for our 
case are the kinematics inside the half-light radius (curves with 
< 50% light/mass enclosed in figure ITTb - The system relaxes 
very quickly in the central regions. After six relaxation times the 
inner regions of the cluster are almost isotropic with ft values 
below 0.2. The relaxation time of NGC 6388 is of order 10 9 yrs 
(at the half-mass radius, lHarrislll996l) which implies that the 
cluster is about 10 relaxation times old. This shows, that high 
anisotropics (such as the the ones needed from our best fit mod- 
els) in NGC 6388 are not stable and would have vanished over a 
short time scale. Therefore, we limit our discussion to the results 
from our isotropic models. 

5.4. Remaining concerns 

In our analysis we used the assumption of a constant M/Ly for 
the entire cluster. In reality, this ratio can vary with radius and 
it may not be well described with a single value if the cluster 
is mass segregated. The rise in the velocity-dispersion profile 
could then be mimicked by a dark remnant cluster at the center. 
Such a scenario is expected in core collapsed clusters which dis- 
play surfa ce brightness profiles w ith logarithmic slopes of -0.8 
or higher (Baumgar dt et all2005l) . We looked at the deprojected 
light-density and mass-density profiles resulting from the sur- 
face brightness profile and the velocity-dispersion profile. The 
mass-density profile from the kinematics drops with oc r~ 2 which 
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Fig. 10. The result of the Monte Carlo simulations of the surface 
brightness profile. For each fixed /3 value 1000 realizations of the 
light profile have been performed. The black dots mark the av- 
erage value of the best fit black-hole mass. The shaded contours 
represent the 68, 90 and 95 percent confidence limits. 
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Fig. 11. The anisotropy /3 as a function of time for different re- 
gions in the cluster. The different colors represent different per- 
centages of the cluster total light/mass going from the inner 
part of less than 10% to the outskirts of the cluster with 90% 
light/mass enclosed. 



would be expected from a core collapsed cluster. But the light 
profile is shallow and does not support that hypotheses: with 
a concentration of c = 1 .7 and a logarithmic slope of the sur- 
face brightness profile of a = -0.28 + 0.08 NGC 6388 does not 
show any signs of core collapse in the distribution of the visible 
stars. We compared these profiles with simulated core collapsed 
or ma ss segregated clusters (described in Baumgardt & Makino 
2003) and were not able to find a good agreement with either 
shape of the light profile or slope of the mass density of the clus- 
ter. We currently consider the presence of a dark remnant cluster 
unlikely but we will perform detailed n-body simulations with a 
variable M/L in the near future. 



6. Summary and conclusions 

We derived the mass of a potential intermediate-mass black hole 
at the center of the globular cluster NGC 6388 by analyzing 
spectroscopic and photometric data. With a set of HST images, 
the photometric center of the cluster was redetermined and the 
result of NG06 confirmed. Furthermore, a color magnitude dia- 
gram as well as a surface brightness profile, built from a com- 
bination of star counts and integrated light, were produced. The 
spectra from the ground-based integral-field unit ARGUS were 
reduced and analyzed in order to create a velocity map and a 
velocity-dispersion profile. In the velocity map, we found sig- 
natures of rotation or at least complex dynamics in the inner 
three arcseconds (0.15 pc) of the cluster. We derived a velocity- 
dispersion profile by summing all spectra into radial bins and 
applying a penalized pixel fitting method. 

Using the surface brightness profile as an input for spheri- 
cal Jeans equations, a model velocity-dispersion profile was ob- 
tained. We ran several isotropic models with different black-hole 
masses and scaled them to the observed data in order to measure 
the mass-to-light ratio. Using y 2 statistics, we were able to find 
the model which fits the observed data best. We run Monte Carlo 
simulations on the inner points of the surface brightness profile 
in order to estimate the errors resulting from the particular choice 
of a light profile. From this, we determined the black-hole mass 
and the mass-to-light ratio as well as the scatter of these two re- 
sults. The final results, with 68% confidence limits are: A black- 
hole mass of ( 1 7 ± 9) x 1 O 3 M and a global mass-to-light ratio of 
M/Ly - (1. 6+0. 3)M /L . In addition, we run anisotropic mod- 
els. The confidence of black hole detection is decreasing rapidly 
with increasing radial anisotropy. However, using N-body simu- 
lations, we found that in a relaxed cluster such as NGC 6388 an 
anisotropy higher than /3 = 0.2 inside the half mass radius will 
be erased within a few relaxation times. 

With respect t o the black-hole mass estimate of 5.7 x 10 3 M 
by lLanzoni et al.l (|2007), which they derived from photometry 
alone, our derived black-hole mass is larger by a factor of three 
with our surface brightness profile profile. Using their light pro- 
file, we obtained a black-hole mass of M. =(11+ 5 ) x 1O 3 M 
which includes their value within one sigma error. ICseh et all 
(2010) obtained deep radio observations of the inner regions 
of NGC 6388 and discussed different accretion scenarios for a 
possible black hole. They found an upper limit for the black- 
hole mass of ~ 1500M o since no radio source was detected at 
the location of any Chandra X-ray s ources. However, inserting 
our black-hole mass in equation 5 of ICseh et al.l (1201 Ol) and us- 
ing their assumptions as well as the Chandra X-ray luminosity 
of the central region, results in minimal but possible accretion 
rates and conversion efficiency (etj ~ 10~ 7 ). Considering the fact 
that supermassive black holes have mass es not much higher than 
0.2 % of the mass o f their host systems (Marco ni & Huntl 2003: 
lHaring & Rixl 12004). our mass with 0.9% seems to be higher 
than the predictions for larger systems. Globular clusters in con- 
trast lose much of their mass during their evolution. This could 
naturally result in higher values of black-hole mass - host system 
mass ratios. Due to the complicated dynamics, the 1 <x uncertain- 
ties are 40 % for our black-hole mass and 10 % for the derived 
M/L v . Figure [T2l shows the position of NGC 6388 in the black- 
hole mass velocity dispersion relation. With our derived effec- 
tive velocity dispersion of cr e = (18.9 + 0.3) kms _I the results 
for NGC 6388 seem to coincide with the prediction made by the 
M. - cr relation. 

The mass-to-light ratio of M/L v - (1.7 ±0.3) M e /L Q derived 
in this work is consistent within the error bars with the dynami- 
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cal results of iMcLaughlin & van der Marell (120051 M/L v ~ 1 .8) 
but slightly lower than the value predicted by stellar population 
models. According also to McLaughlin & van der Marel (2005), 
we would expect a value of M/L v - (2 .55 + 0.28) Mq/Lq for 
[Fe/H] = -0.6 dex and (13 ± 2)Gyr. iBaumgardt & Makinol 
(2003) have shown that dynamical evolution of star clusters 
causes a depletion of low-mass stars from the cluster. Since 
these contribute little light, the M/Ly value drops as the clus- 
ter evolves. Dynamical evolution could therefore explain part of 
the discrepancy between theoretical and observed M/Ly values. 

Summarizing, it can be said that the globular cluster NGC 
6388 shows a variety of interesting features in its photomet- 
ric properties (e.g. the extended horizontal branch) as well as 
in its kinematic properties (e.g. the high central velocity dis- 
persion or the rotation like signature in the velocity map). In 
this work we investigated the possibility of the existence of an 
intermediate-mass black hole at the center. With simple ana- 
lytical models we were already able to reproduce the shape of 
the observed velocity-dispersion profile very well if the clus- 
ter hosts an IMBH. Future work would be to compute detailed 
N-body simulations of NGC 6388 in order to verify whether a 
dark cluster of remnants at the center would be an alternative. 
Additionally, it is crucial to obtain kinematic data at larger radii 
to constrain the mass-to-light ratio and the models further out. 
Also proper motions for the central regions would further con- 
strain the black hole hypothesis. 

The study of black holes in globular clusters is currently lim- 
ited to a handful of studies. For this reason, it is necessary to 
probe a large sample of globular clusters for central massive ob- 
jects in order to move the field of IMBHs from "whether" to "un- 
der which circumstances do" globular clusters host them. Thus 
tying this field to the one of nuclear clusters and super-massive 
black holes. From the experience we gathered with NGC 6388, 
we can say that the dynamics of globular clusters is not as simple 
as commonly assumed and by searching for intermediate-mass 
black holes, we will also be able to get a deeper insight into the 
dynamics of globular clusters in general. 
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